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Abstract 

A new theory of hydrodynamics of uniaxial nematic liquid crystal films 
in the presence of defects is developed. A gauge field incorporating screen- 
ing is introduced, resulting in the static elastic free energy having the form 
of a two-dimensional Abelian-Higgs model. Hydrodynamic equations are 
derived via the standard methods of de Groot and Mazur. By working in 
the vicinity of the Bogomol'nyi equations consequences for defect centre 
motion are outlined. 

PACS numbers: 61.30.Jf, 61.30.Cz, 11.15.Kc 

Liquid crystals are important for technical applications such as liquid crystal 
displays. However the dynamics of defects studied in experiments is also used 
to model the evolution of defects like strings in cosmology [jij . Here we propose 
a theory to model the dynamics of nematic liquid crystals with defects. 

1 Description of Defects 

Liquid crystals are states of matter which are in between those of a liquid and 
crystal. In terms of their mechanical properties they are akin to isotropic fluids, 
but in terms of their optical properties they behave more like crystals. There 
are various types of liquid crystals. We shall only consider uniaxial nematics. 
The nematic molecule is a rod-like molecule and the average direction of these 
rods is given by a vector n with fixed length, n = 1. The distribution of the 
molecules is random, but there is orientational order. For this macroscopic 
continuum theory the molecules are like head- less arrows, so n and -n have to 
be identified since they represent the same physical state |2|, Q . Here we study 
nematic liquid crystals in two-dimensional systems. These are thin layers, or 



'Based on lecture by S.S. at XXXVIII Cracow School of Theoretical Physics, Zakopane, 
Poland, June 1-10, 1998. 

^E-mail: gk@maxwell.ph.kcl.ac.uk, ss2@maxwell.ph.kcl.ac.uk 



1 



Abelian-Higgs-Navier-Stokes Hydrodynamics ... — G. Kurz, S. Sarkar 2 

films, where the molecules align parallel with the surface layer, i.e. the director 
is confined to two dimensions. 

The symmetry of possible rotations of the director is spontaneously broken, 
since the molecules tend to align themselves in parallel with one another. This is 
reflected in the form of the distortion free-energy density describing the statics 
of nematics. A ground state which minimizes this density is locally aligned 
(with possible exceptions at singular points or lines, also called defects) @. 
The most general form of the distortion free-energy density is 

4> d = \K X [V ■ n] 2 + lK 2 [n ■ (V x n)] 2 + ±K 3 [n x (V x n)] 2 . (1) 

The positive constants K%, K 2 , and K% are associated respectively with splay, 
twist, and bend deformation. For a classic nematic like p-azoxyanisol (PAA) at 
120°C their values are K\ = 0.7 x 10 6 dyn, K 2 = 0.43 x 10 6 dyn, K 3 = 1.7 x 10 6 
dyn Q. As (f)d is even in n, the states n and -n are indeed indistinguishable. 
We remark that in two dimensions the three kinds of deformation are not in- 
dependent, since there we have 

[n • (V x n)] 2 = [V • n] 2 + [n x (V x n)] 2 . 

For our purpose we use a simplified version of (ftd with enhanced symmetry, the 
one-constant approximation where one assumes K\ = K 2 = K% =: K 

4>i = \Kdjnidjni . (2) 

(The convention of summing over two repeated indices is used throughout, 
unless mentioned otherwise.) 

The ground state of the nematic is described by a direction or orientation. 
Consequently in the three dimensions the manifold of such states (known also 
as the vacuum manifold) has the topology of a sphere S 2 . However, since n and 
— n are identified, this manifold is further restricted to the projective plane. 
In two dimensions the manifold is a circle with opposite points on a diameter 
identified, which is still homeomorphic to a circle, S\ 

The possible topological defects are given by the homotopy groups of the 
vacuum manifold. For S\ the only non-trivial homotopy group is tti(Si) which 
determines the possible types of defect points (also known as disclinations). 
These points can be considered as strings, if we see our two-dimensional sys- 
tem as the cross-section of a three-dimensional one. Other homotopy groups 
for systems with different vacuum manifolds determine walls, monopoles and 
textures; a brief introduction is given in tti(Si) reflects the number of topo- 
logically distinct noncontractible loops on Si. This homotopy group can be 
represented through integers, or more precisely through half- integers. Defect 
points are described mathematically like vortices, and in nematics they can 
have half-integer winding numbers because of the identification of n and — n. 
Consequently the defect points in two-dimensional nematic liquid crystals are 
characterized by their winding numbers =L^, ±1, ±|, . . . Schematic illustrations 
of director configurations and pictures from experiments, both for various types 
of defects, can be found e.g. in f| §, §, g. 
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To deal with the possible singularities of the director, it has been found 
useful to soften the constraint n = 1. The length of the director is regulated 
instead by a potential term which is added to the energy density. We take the 
potential to be proportional to (|n| 2 — l) 2 . This is minimized for |n| = 1. So 
the ground states remain the same, and so does the vacuum manifold of the 
system. 

In the past the introduction of gauge fields has been suggested and used to 
model systems with defects |7|, ||, ^, [jl], |llj . A Gauge field can model screening 
when many vortices are present. Here we propose a new form of the distor- 
tion free-energy density which will be used later on to derive the equations of 
hydrodynamics, in the same way as (ftd was the basis for the Ericksen-Leslie 
equations which describe the hydrodynamics of nematics. It is recognized that 
the Ericksen-Leslie equations do not satisfactorily incorporate disclinations |||. 

For two-dimensional liquid crystals the relevant gauge theory with the vac- 
uum manifold Si is the Abelian-Higgs model. This is a relativistic model for- 
mulated for three space plus one time dimension. We take on just the static on 
space part. It gives us a generalized form of the one-constant approximation 
01 in (0). When written in terms of the complex director n = n% + in-i relating 
to the two-dimensional system, it has the form 

<j> = -KDjnljjn + \X(\n\ 2 - l) 2 + \F 2 . (3) 

The ordinary derivatives in cj)\ are replaced with covariant ones, Dj = dj — igAj, 
and a term in the derivatives of the gauge field A, the gauge field strength 
F = d\A<i — 82 A\, is added. We also have the potential mentioned above. 
Kawasaki and Brand [|ll]] claim to derive a free-energy density close to (||) 
from (||) by carefully separating the director into singular and non-singular 
pieces. These arguments are delicate especially those involving the number of 
the degrees of freedom of the system. 

Whereas 0i is invariant under global rotations %(r) i— > Rjkn>k(r), with Rjk 
being a constant rotation matrix, <j) is covariant, i.e. invariant under a local 
rotation of the form 



n(r) i — y e iv(r) n(r) , 
Mr) ~ AjW + i^r). (4) 

This is referred to as a gauge transformation determined by the space-dependent 
phase angle <^(r). 

The covariant energy density (||) implies that the interaction between discli- 
nations is screened. Mathematically the decay of the force between two defects 
in a model without gauge fields is inversely proportional to the distance be- 



tween defects [12]. In the Abelian-Higgs model the force falls off exponentially 



with the distance 13, 14 1 so it has a weaker long range interaction. Even 
though very little is known about the dynamics of disclinations in nematics § 
it is still expected that a large number of them causes screening and that they 
influence dissipative flow properties. Our model proposes to take account of 
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this by introducing the gauge field, which has also topological and geometrical 
significance. 

2 Hydrodynamic equations for Nematics with de- 
fects 

We derive hydrodynamic equations to model nematics with defects based on the 
covariant energy density (^). We will follow the standard procedure of hydro- 
dynamics based on the concept of the fluxes and forces of de Groot and Mazur 
p~5[] . This procedure was successful in deriving the usual dynamics of nematics. 
The resulting equations are referred to as the Ericksen-Leslie equations and the 
variables are the velocity field v and the director n. On suppressing n they 
reduce to the Navier-Stokes equations. In our hydrodynamics for defects in 
addition to v and n, the gauge field A is introduced. Consequently not only 
additional terms are introduced into the existing equations, but there is a new 
dissipative equation for A. 

In this section we will employ the notation n = (711,712) and A = (A\,A2) 
for the two-dimensional vector fields. The components are functions of the 
two-dimensional coordinate r. Instead of the complex number i the covariant 
derivative contains the matrix I which is given by 




(I is the generator for the group of rotations in the plane). The covariant 
derivative is then Dj ^i = djS k i — gAjI k [. Its action on n we shall also note 
in short as (Dn)j k = Dj^rii. This expression translates into components the 
complex one used in equation @ by making the real and imaginary part the 
first and second component in the second index in {Dn)j k . Furthermore, the 
expression for the gauge field strength Fj k = djA k — d k Aj is used instead of F. 
Our theory is dissipative and does not have the relativistic invariance of field 
theories. Consequently time and space components are not on an equal footing, 
and we are not forced to introduce an Aq field. We will consider Aq = and 
time independent gauge transformations. The covariant form of the distortion 
free energy density for nematics from (^) can then be rewritten as 

= \K{Dn) jk {Dn) jk + \\{n jnj - if + \F jk F jk . (5) 

Physically the system is understood as a layer with a small but homogeneous 
thickness, and the vector fields are confined to the plane, as indicated above. 
The energy of the whole system is written as <3? = J d 3 x 0. 

Instead of equation (|j) the gauge transformation, under which the energy 
density is invariant, now looks like 

rij(r) i-> R jk (r)n k {r) , 
A^Im ^ AjWlM + ^djRkiWR-^r) , (6) 
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where Rjp.(r) is a two-dimensional rotation matrix. As a consequence Fji~ is 
again invariant and the covariant derivative of the director transforms like 
(Dn) jk (r) i-> Rki{r) (Dn)ji(r). 

The dynamics of nematics requires the following two variables appropriate 
for isotropic fluids in addition to the director n(r) and the gauge field A(r): 
the fluid velocity v(r) = (vi(r),V2(r)) and the pressure p(r) (discussed further 
below). The density p of the fluid is assumed to be constant. This implies the 
incompressibility condition 

d jVj = . (7) 

The acceleration of the fluid is given by the derivative of the tensor of the total 
stress <x(r) (by Newton's law combined with the definition of the stress tensor 
in dfi = d 2 Sj (jji, which the force dfi on a surface element d 2 Sj), 

pd t Vj = dicrij . (8) 

The total stress is split into the elastic stress cr e (also called the Ericksen stress) 
and the viscous stress tr v , 

a = a e + a" . (9) 

In addition to ([?]) and (||) the dynamics of nematics requires equations for the 
first-order time derivatives of the director and gauge field. These, together with 
the expressions for elastic and viscous stress as functions of the above mentioned 
variables, are derived in what follows. 

2.1 The Entropy Source 

From the expression of the entropy source for a flowing nematic liquid the hy- 
drodynamic fluxes and forces are identified. For small fluxes the forces are linear 
in the fluxes. The aim is to find these linear equations or constitutive relations. 
Their proportionality factors are of the form of a viscosity coefficient multiplied 
by a monomial in the hydrodynamic variables. The viscosity coefficients need 
to be determined phenomenologically. 

The dissipation (product of temperature and entropy source) is the decrease 
in stored free energy. Its contributions are the kinetic and distortion energy. 
Any external fields (electric, magnetic or gravity) are not taken into account. 

TS = -j t Jd 3 x(±pv j v j + <t>) , (10) 

where S is the entropy. 

To calculate the first term of the right hand side we use the fluid acceleration 
equation @. On using integration by parts we obtain 



d 
dt 



J d 3 xpvjVj = J d 3 x (J jk djVfc -\- surface terms . (11) 

To calculate the second term of the right hand side in (jlOj ) we perform 
two types of variations of the distortion energy density. One is a variation of 
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the energy density by changing the orientation of the director and the gauge 
field. The transformations are n i— > n' + 5n and Ah A' + 5 A. Through partial 
integration we get for the energy 

5$ = J d 3 x (-h] 5nj - hf 5A^j + surface terms (12) 

with 

h " ~ B 'wk) ~ 3& " K(DDn) <" ~ A( "<"< " ■ <13> 

= *5(^) " a£ = d,F " - Ks(Dn) ' 1 '» "' ■ (14) 

Usually, h n is called the molecular field; here we mark it as related to the 
director n. The additional expression which appears in our model has the 
analogous notation h A and is called the molecular field related to the gauge 
field A. The double application of the covariant derivative in (|l3| ) is explicitly 
given by (DDn)i m j = Dij^D^^m. We notice that the two molecular field are 
non-zero only out of equilibrium, as the Euler-Lagrange equations belonging to 
(f> are h n = and h A = 0. 

The second variation of the energy density consists of displacing the liquid 
without changing the director orientation. The variation parameter is u(r) 
which gives the displacement of the coordinate r to r' = r + u(r). The director 
transforms such that 

n(r) h- > n (r ) = n (r + u) = n(r) . (15) 

The transformation of r to r' causes the derivative operator to transform as 
dj i — > d'j = (Sjk — djUk(r))dk- Gauge invariance is obtained by imposing the 
same transformation on the gauge field. All calculations are just up to linear 
order in u. This gives the following results for the variation of the variables: 



6 (djUk) = dju'^r') - djn k (r) = -din k djUi , 

5A 3 = A' 3 {v') - Aj(r) = -A, S,-«, , (16) 
5 (djA k ) = d'jA'^v') - djA k (r) = -d t A k d jUl - d J A l d kUl , 

and by using the anti-symmetry of Fj k in the term involving 5(djA k ), we get 
the variation of the energy. The result is written in terms of the elastic stress 
cr e , which is the usual variable obtained from this kind of variation, 

5$ = J d 3 xa e jk djU k , (17) 

a% = —K(Dn)ji(Dn)hi - F 3l F kl - P 5 jk . (18) 

The function p = p(r) is a Lagrange multiplier that enforces the incompress- 
ibility condition (0) on the displacement given by u(r). It is called pressure 
and has the negative sign reflecting that generally the stress tensor a for an 
isotropic fluid without friction reduces to aj k = —Sj k p where p is the ordinary 
pressure. 
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The total variation is obtained by adding ( |l2|) and (17). This gives the total 
time derivative 



4: J d 3 x cp = J d 3 x (— c| fc djV k + fij + hf Aj^j + surface terms , (19) 



dt 



where u was replaced by v, ri is understood as the material time derivative, 
the rate of change of the director as experienced by a moving molecule, and is 
written in terms of covariant quantities, 

hj = (Dn) 0j +vi(Dn) lj . (20) 
Analogously A is defined by 

A j = F 0j + vi . (21) 

Under gauge transformations n behaves like n, and A is invariant. 

Finally we get the entropy source by putting ([ll]) and ( p~0| ) into equation 

©■ 

TS = J d 3 x crj fc djV k + h" hj + hf Aj + surface terms , (22) 

where equation (§) was used to replace the total and elastic stress by the viscous 
stress. 

2.2 The Balance of Torques 

From the entropy source (^) the hydrodynamic fluxes and forces could be 
identified already. However we continue to follow de Gennes' derivation and 
rewrite it further by splitting the viscous stress into a symmetric and an anti- 
symmetric part. 

°jk = \ (vjk + Vjk) , ( 23 ) 

r i = - e m a ii ' ( 24 ) 

where ejki is Levi-Civita's anti-symmetric symbol. The same is done for the 
tensor formed by the derivative of the velocity, 



s 



ik = I (djV k + d k Vj) , (25) 
w i = l e jki dkVi ■ (26) 

(For the cross products the variables like v and cr v are seen as three-dimensional 
where the components relating to the third dimension are zero. As a result T 
and u are perpendicular to the plane of our two-dimensional system.) This 
allows us to write 

a v jk djV k = crj| s jk - Tj uj . (27) 
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We will now show that T may be identified as the torque density and can 
be written as T = n x h n . To further this aim the first step is a variation of 
the energy density for a constant director rotation around the axis parallel to 
5u with a small rotation angle \8u)\, 

n i— > n + 5u x n . (28) 

This is a gauge rotation as given in equation (||), but with a constant rotation 
matrix, and so the corresponding variation of the gauge field is zero. By sepa- 
rating into volume and surface terms and using the molecular field h n of (|b|), 
we get a change of the energy <I>y = f v d 3 x (ft of any part V of our system given 
by 



5$ 



v 



/ d 3 x W e jM m+ d 2 Si 
Jv J Jav 



d(d i n j 



Su k (29) 



due to the variation. The second integral inside the bracket is the torque on the 
director on the surface. Since the energy density is invariant under the gauge 
rotation (|28|), we have 5<&y = 0, and the torque on the director can be replaced 
by the first integral inside the bracket. The second contribution to the torque 
comes from the stress on the surface through J ej k i r k dfi = J ej k i r k (Tii d 2 Si, and 
so the total torque on the volume V is 

Tj = / d 3 xejkin k hf + / e jk ir k au d 2 Si . (30) 
Jv Jav 

On the other hand the torque is given by the rate of change of the angular 

momentum L, i.e. L = T. It is calculated from L = J v d 3 xpr x v by putting 

in the total stress from @ , integrating by parts and by using the equations (||) 

and (^) to rewrite the volume integral. We then get 

Lj = / e jk ir k aud 2 Si+ / d 3 x (Tj - e jkl a kl ) . (31) 
JdV Jv 

The term in <r e is zero. This is based on the fact that the energy density is invari- 
ant, i.e. 5&v = 0, under coordinate rotation r i— > r + 5u x r (without changing 
the director) around the axis 5u with a small rotation angle. This variation is 
the same as that carried out in equations ( |l5|) to ( |l8| ) with u(r) = u> x r. (It is 
also valid for any part V of the system.) Hence <5<3>y can be obtained from (|D 
where djU k = e jk ioJi. 



Comparing equations (30) and ( |3l|) and using the fact that they are valid 



for any volume V within the system, we can conclude that 



r = n x h n . (32) 



By using equation (|27| ) and (^) in (|22j ) , the new form of the entropy source 
is obtained as 

TS = J d 3 x af k s jk + h n j Nj + hf Aj + surface terms , (33) 
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where we have introduced the vector N, the rate of change of the director with 
respect to the background fluid. 



N = n- w xn. (34) 

Since u is perpendicular to the plane of the film containing n, uxn behaves 
in the same way as n under gauge rotation (and so does ri); so N rotates like 
n for a change of gauge. 

The entropy source ( |33|) displays three types of dissipation: dissipation by 
shear flow, by rotation of the director with respect to the background fluid, and 
by change of the gauge field. 



2.3 The Linear Equations between Fluxes and Forces 



In the entropy source (|33|) each contribution is identified as a product of a flux 
with the conjugate force. We choose them in the same way as de Gennes Q: 
sn, si2, N, and A as hydrodynamic fluxes conjugate to the forces of/, 2cr^|, 
h n , and h A respectively (the factor 2 reflects that of! s±2 occurs twice). In the 
limit of weak fluxes, the forces are linear functions of the fluxes. We can express 
this in the form 



V 



9rr vs 
zo 12 

hi 



\ 



( 



) 



\ 



constituent 
matrix 



V 



/ 



N p 

V 4 J 



(35) 



A priori this matrix is a tensor constructed from the variables of the prob- 
lem. This gives the matrix a large number of possible terms. However, certain 
conditions restrict the form of the matrix. We look at the behaviour of the 
fluxes and forces under time reversal. From equation (|8|) and ( |i~8|) it can be 
seen that the components of cr and cr e are even under time reversal, and so 
er v is even too. The same goes for the other hydrodynamic forces, whereas all 
the fluxes are odd. The Onsager reciprocal relations [15] then imply that the 
constituent matrix is symmetric. We demand that the matrix be compatible 
with the local symmetry of uniaxial nematics, Dqo^. This excludes derivatives, 
the gauge field and certain tensor types from the matrix. Furthermore the ma- 
trix elements need to be such that the equations are invariant under director 
inversion n i— ► — n, which is essential for nematics. (This inversion leaves the 
gauge field unchanged, since n i— > — n corresponds to a gauge transformation 
with constant phase angle tp = 7r.) The only scalar invariant is |n| 2 , but outside 
the core of a disclination it is of value one and near the defect centre it is small; 
like in the Ericksen-Leslie theory it will not be used in the constituent matrix. 

With these conditions ( |35|) translates into the following equations 



cr?? = ai TiiUjUkniSki + ^72 (riiNj + rijNi) (36) 

+ «4 Sij + \ (a 5 + a 6 ) (riiSkj + rijSki) n k , 
hV; = 7iiVi + ^njSji , (37) 
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hf = /Mi + P 2 n i n j A j . (38) 

The viscosity coefficients oci in cr vs are set up with the result for o~ v in mind. 
From equations (23) and ( |24|) we know that of- = <r? ? — \eijkTk- Hence with 



the help of equations (|36l ) , (|32| ) , and (^) we obtain the total viscous stress 



ai niUjUkniSki + a 2 niNj + a 3 rijNi (39) 

+ «4 Sij + Q 5 UiUkSkj + 0*6 njHkSki , 



where a 2 and 03 are defined such that 71 = 03 — a 2 and 72 = «3 + «2 • (In 
addition we have the relation 72 = 01% — a§ j2j.) The coefficients on and % 
are viscosities known from the Ericksen-Leslie theory M. The coefficients Pi, 
related to the dissipation of the gauge field, are new and have been introduced 
by us. 

In summary we have the hydrodynamics of nematic liquid crystals with 
defects given by the three dissipative equations (||), fl37|), (p8|), corresponding 
to the three hydrodynamic variables, the velocity field v, the director field 
n, and the gauge field A. In addition we have the incompressibility condition 
(0) giving rise to the pressure variable p. To specify the three hydrodynamic 
equations in detail we need the total stress a in @, the elastic stress er e in 
(|l8|), the viscous stress cr v in (|39"|), the rate of change of the director N in 



the material time derivatives ri and A in (|2C]) and ( 21 ) , and the molecular fields 
h n and h A in @ and (0). 

If the gauge field is zero, the original form of the Ericksen-Leslie equations is 
obtained for the director and velocity field from (||) and (p7|), together with the 
incompressibility condition. Since the original equations due to the dissipation 
were not invariant under a global rotation, the new equations are not invariant 
under gauge transformation. If we require gauge invariance we can keep only 



the diagonal terms in the constituent matrix in (^). Instead of (37) and ([" 
the dissipative equations for the director and gauge field become 



K = liNi, (40) 
hf = frA, (41) 

and the viscous stress from (^) reduces to 

Oij = 571 faNi - riiNj) + a 4 . (42) 

(The 71 term is invariant under gauge rotation because we work in only two 
dimensions.) 

In a first attempt to solve these equations we truncated the system by 
assuming zero velocity, i.e. the configuration of the system changes through re- 
orientation of the molecules without displacing them. We then have to consider 



just the two dissipative equations (|40[) and (|4l|), and this is a simplification 
used earlier flTBl. 
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3 Hydrodynamics of defects in the Abelian-Higgs 
model 

In this section we summarize very briefly how the general equations for the 
hydrodynamics of nematics with defects derived in the previous section are 
used to find the dynamics of disclination points in a simplified system with zero 
velocity field. A detailed account of the related calculations and results can be 



found in [17 1 . 



First we consider the distortion free-energy density (y). For a certain value 
of the gauge coupling g, the Bogomol'nyi limit, a configuration of disclinations 
with winding numbers of the same sign is a static solution. In general, of course, 
several defects are not a static solution, and there is a force between them which 
makes them move. However if the coupling deviates only slightly from the Bo- 
gomol'nyi limit we can still assume that the configuration takes on the form 
of a static solution at every instant. Such a configuration change over time is 
known as quasi-static motion and is induced by the force due to the deviation 
from the Bogomol'nyi limit. A solution with defects in the Bogomol'nyi limit 



can be parameterized by the defect positions [18|. The motion of the config- 
uration is then given equivalently by the motion of the defect positions in the 
space of possible static solutions (also known as moduli space). The concept 



of quasi-static motion was developed by Manton |19[ and has been applied by 
Samols and Dziarmaga |2l], p2| . 

Physically a defect is a macroscopic region with a certain topological charge. 
Mathematically this is modelled as a topological defect point with a core region 
and a well defined centre point (providing the defect position coordinate). If 
the distance between the defects is large compared to the core size, a static 
multi-defect solution in the Bogomol'nyi limit is found to a good approximation 
by superposition of solutions with single defects in different positions. A static 
solution for several defects with overlapping core regions is modelled by a single 
defect solution which is modified such that it allows for disintegration of this 
defect into smaller ones. 

In both cases, for large distances or overlapping cores, the static solution 
in the Bogomol'nyi limit can be parameterized by the defect positions. Their 
dynamics is given by the field equations (^) and (41). These are reduced 
to equations just for the defect positions by projection on the moduli space. 
The moduli space is spanned by zero-modes of the linear operator given by 
linearization of equations ( f40|) and ( fH|) in the field variables n and A. It turns 
out that the zero-modes are derivatives of the defect solutions with respect to 
the defect coordinates. The result of the projection is a set of equations in the 
defect positions and their first time derivatives. 

We derive and solve these equations for both cases mentioned above. For 
defects with overlapping cores they indicate that an unstable disclination with 
higher winding number can disintegrate into smaller ones which move away from 
one another radially with exponentially increasing speed (as long as their core 
regions overlap). Two disclination far apart from one another move on a straight 
line, where their distance increases logarithmically with time. Disclinations with 
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smaller winding number have faster motion. 
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